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A survey of topics of recent interest in Hamiltonian and Lagrangian dynamical 
systems, including accessible discussions of regularization of the central force prob- 
lem; inequivalent Lagrangians and Hamiltonians; constants of central force motion; a 
general discussion of higher-order Lagrangians and Hamiltonians with examples from 
Bohmian quantum mechanics, the Korteweg-de Vries equation and the logistic equa- 
tion; gauge theories of Newtonian mechanics; classical spin, Grassmann numbers, 
and pseudomechanics. 

I. INTRODUCTION 

The study of classical mechanics is vast and ancient. Therefore, this collection of results, 
observations and questions necessarily omits most of the field and probably misses a num- 
ber of older references even for topics covered in detail. We focus principally on issues of 
symmetry and subjects (old and new) which have appeared in the literature within recent 
decades. Nor should it be thought that we provide a complete survey of even the results we 
do discuss. Instead, our references for each topic are probably sufficient for the interested 
reader to gain a foothold on the relevant research. 

Of course, the list of topics we do not examine is extensive. Certain topics such as 
nonlinear dynamics receive only a brief mention as an example of a higher order system in 
Section 5. We have chosen to omit any discussion of electromagnetism while touching on 
special and general relativity only to provide examples. Since our presentation is intended 
for a broad audience, we have also avoided the large body of formal work. Thus, while the 
formal study of symplectic manifolds, Kahler manifolds, Poincare sections and so on make 
heavy use of modern differential geometry and field theory techniques, little mention is made 
of progress in these directions. 

What remains is nonetheless filled with fascinating and diverse surprises in a field often 
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mistaken to be complete. Thus, what we do cover is a wide array of topics ranging from the 
Kepler problem to supersymmetry. The unifying theme, if there is one, is the occurrence 
of unexpected and surprising symmetries in classical physics, and especially in Lagrangian 
and Hamiltonian dynamics. Though we treat a few topics simply because there is recent 
reference to them in the literature, most of the topics concern symmetry in one way or 
another. The uses vary greatly, from the use of anticommuting numbers to the rotation 
group, from unusual constants of motion of the Kepler problem to the infinite hierarchy 
of constants of motion of the Korteweg-de Vries (KdV) equation. An additional guiding 
principle has been to treat topics that may not be familiar to many readers. 

Curiously, quantum mechanics and quantum field theory have had a strong impact on cur- 
rent work in classical physics. As a result, a brief discussion of quantum mechanics appears 
in our examination of inequivalent Lagrangians in Section 3 and Bohmian quantum mechan- 
ics is discussed in Section 5. Further connections between quantum and classical mechanics 
are suggested in our treatment of Lagrangian and Hamiltonian dynamics as gauge theories 
in Section 6. Finally, Section 7 owes its entire existence to insights from supersymmetric 
quantum field theory. 

The organization of the paper is simultaneously from old to new and from easy to diffi- 
cult. Section 2 provides a warm-up exercise with some new thoughts on an old topic - the 
regularization of the central force problem. From there we move gradually to more recent 
and more mathematically challenging questions. In Section 3 we discuss inequivalent La- 
grangians and Hamiltonians, a topic which begins with Lie and Dirac (if not earlier) and 
which received considerable new input in the 70s. Through the same period, old knowledge 
resurfaced with the rediscovery of the Laplace- Rungc-Lenz and Hamilton vectors. A deriva- 
tion of these rediscovered constants of the motion is given in Section 4, using a technique 
based on an old theorem. While the theorem will no doubt be familiar to mathematicians, 
its simple application to finding constants of the motion does not appear in classical physics 
textbooks. 

Moving to more active areas of current interest, we look at the occurrence of higher 
order differential equations in classical physics. After a brief general introduction to these 
systems at the beginning of Section 5 is an example of such equations - Bohmian quantum 
mechanics. 

The final two Sections deal with truly contemporary insights. In Section 6, a development 
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of both Lagrangian and Hamiltonian dynamics as gauge theories shows an interesting new 
connection between classical physics and conformal symmetry. Then, in Section 7, two 
further developments of field theory - spinors and anticommuting variables - are discussed 
in the context of particle mechanics. The KdV equation and the approach to chaos are 
treated in Appendices. 

Before embarking, a few general comments are in order. First, observe that each section 
below is essentially independent of the others. Each section has its own brief introduction 
and references. Note that the end of most sections we have tried to provide a few stimulating 
questions. These questions do not reflect any consensus thinking and should not be taken 
to be the definitive puzzles facing the field. Rather, they are suggestions of some directions 
which might or might not prove fruitful. Finally, it should be noted that where derivations 
are given without citation, we have produced original calculations. We make no further 
mention of this fact since it is probable that many or all of these calculations have already 
appeared somewhere within the last few hundred years! 



II. REGULARIZATION OF THE CENTRAL FORCE PROBLEM 

We begin with some of the oldest problems of classical physics. In this Section and the 
next, we explore some interesting features of the Kepler problem and other central force 
motion. In this Section, we examine regularizations of the Kepler problem. In the next 
section we present a technique for finding constants of the motion then, as an example, 
use the technique to find some recently rediscovered constants of the Kepler problem ^. 

Regularizations of dynamical problems are transformations that turn the equations of 
motion into a simpler or less singular mechanical problem. Euler Q] and Levi-Civita [4J pro- 
duced one- and two-dimensional regularizations, respectively, of the Kepler problem. These 
regularizations turn the Kepler /Coulomb equation of motion into an isotropic oscillator. It 
is not surprising that this is possible, because the transformations are time-dependent. In- 
deed, using similar transformations, it is possible locally to turn any central force problem 
into the isotropic oscillator. We present a general proof of this claim below. The Euler and 
Levi-Civita results are special cases. 

There are some recent discussions in the literature extending these regularizations. Since 
the Levi-Civita result makes use of a complex variable, some authors have explored the 
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idea that the use of a vector space which is a. 



so a number field gives insights into the 



problem. Thus, Kustaanheimo and Stiefel ([^, give a quaternionic transformation 

from the 3-dim Kepler problem to a constrained 4-dim isotropic oscillator, therebv show- 
ing that bounded Kepler orbits have an underlying SO (A) symmetry. Bartsch js| writes 
the Kustaanheimo-Steifel result in terms of Hestenes' geometric algebra 0]. Such use of 
quaternionic, Clifford or Grassmann variables (see below) often extends, or streamlines, the 
presentation of classical results. 

However, it seems unlikely that number fields are necessary to transform the Kepler 
problem into the oscillator. If that were the case, we would expect regularization to be 
possible only using real, complex, quaternionic or octonionic variables and therefore only to 
occur in dimensions less than or equal to eight. But since both the Kepler problem and the 
isotropic oscillator are inherently two dimensional, the Levi-Civita solution should suffice in 
any higher dimension as well. Our generalized solution below demonstrates this to be the 
case. 



A. Higher dimensions 



Consider the general central force motion in any dimension d > 2. We begin from the 
action 

S = I dt 



-m— ; V r 

2 dt dt 



where r = ^/Xjxl. It follows that 



dt^ 



-V' — 



We first compute the total angular momentum 



This is conserved, since 



d _ d , 

dt dt^ ^ ^ ' 

( d'^Xk d'^Xj 

= 0. 

To prove from this that the motion hes in a plane, let xq and Vq be the initial position 
and velocity. Then the angular momentum is 

Mij = xoiVoj - xojVoi ^ 

Let w^"^ a = l,...,n — 2, bea collection of vectors perpendicular to the initial plane 

P = {v = axo + /3vo |Vq!, /3} 
w("V = 

so that the set {xq, vq, w*^")} forms a basis. Then, for all a 

wf^M,, = 0. 

Now, at any time t, is given by 

Mij — m {xiVj — XjVi) 

and since M^j is constant we still have 

= wf^^'m (xiVj — XjVi) 

= (w(") • x) V - X (w" • v) 

Suppose, for some ao, that 

^(ao) . X ^ 

Then 

W"o ■ X^ 

and Mjj is identically zero, in contradiction to its constancy. Therefore, we conclude 



V = X 



w 



(«) . X = 
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for all a. A parallel argument shows that 

w^'^) ■ V = 

for all a, so the motion continues to lie in the original plane. 

Now we choose polar coordinates in the plane of motion, and the problem reduces to 
two dimensions. We next need to deal with the presence of angular momentum. With 
coordinates x'-"-' in the w*^"^ directions, the central force equations of motion are 



d r dip dip 
^ ' 'd^~^~dl~dt 



-V'(r) 



d {mr'^ip) ^ 
Jt ~ 

We choose x^""^ = 0, and set L = mr'^ip = constant. Eliminating ip, these reduce to the single 
equation 

d'^r 

m^-^,=-V'ir) (1) 
dt^ mr-^ 

Notice that now any transform of r will change the required form of the angular mo- 
mentum term. What works to avoid this is to recombine the angular momentum and force 
terms. We again start with 

Id 



Then eq.(|T)) becomes 



Rearranging, we have 



f'dr\f'^dT) m^p- df^^^""^^- 



d^u M^r ^,dv^^^ 



dr'^ m^P df 
M^f df dV 



m^P du df ^""^^ 
m?p du du ^"^^ 

To obtain the isotropic harmonic oscillator we require the combined angular momentum and 
force terms to give the needed expression: 

M2 df dV M2 



rn^f^ du du m?u^ 
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Integrating, 



If we define 



(2) 



the required function / is 



-1 I 1 , 2 C 



Substituting this solution into the equation of motion, we obtain the equation for the 
isotropic oscillator, 

3 = 

Therefore, every central force problem is locally equivalent to the isotropic harmonic oscil- 
lator. Of course, the same result follows from Hamilton- Jacobi theory, since every pair of 
classical systems with the same number of degrees of freedom are related by some time- 
dependent canonical transformation. 

The solution takes a particularly simple form for the Kepler problem, V = —a/r. In this 
case, eq.Q becomes 



M2 a f 1 



-ku' 



2m2/2 / \ 2m%2 ' 2 ^" '2 
Solving the quadratic for 1//, we take the positive solution 
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a + 



a2 + 



M2 / M2 



m? \ m?v? 



+ ku^ + c 



1 + 



M 



'ku*+ c + 



amu 



a'^m'^\ „ M2 
H ^ 



M2 



There is also a negative solution. 

We may choose c to complete the square under the radical and thereby simplify the 
solution. Setting 



2\fkMM 



m 



M2 



the positive solution for / reduces to 



1 am^ 
7= 



+ m\/ku+ (M/M)- 
u 
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or 

/= ^ 

{mVk)u^ + {am?/M'^)u + M/M 

The zeros of the denominator never occur for positive m, so the transformation / is regular 
in the Kepler case. The regularity of the Kepler case is not typical - it is easy to see that 
the solution for / may have many branches. The singular points of the transformation in 
these cases should give information about the numbers of extrema of the orbits, the stability 
of orbits, and other global properties. The present calculation may provide a useful tool 
for studying these global properties in detail. The problem of global properties of orbits 
remains open - power law forces have been examined 



but more complicated potentials 



allow arbitrarily many extrema. For example, the potential 



gives an effective potential 



V = a{r — ro)^^ 



Straightforward perturbation about circular orbits shows that, for arbitrary fixed angular 
momentum M, the frequency of radial oscillations may be increased without bound by 
increasing p. Such closed orbits will have arbitrarily many extrema. 



B. Euler's regularization 

Essential features of the regularizing transformation are evident even in the 1-dim case. 
The Euler solution uses the substitutions 

_2 d n d 
x = -u , -7:=u — 
dt dr 

to turn the 1-dim Kepler equation of motion into the 1-dim harmonic oscillator. Before 
moving to a proof for the general n-dim case, we note that more general transformations 
are possible in the 1-dim case. Beginning with the equation of motion, 

d'^x a 
dt"^ x'^ 

let 

„ , , d Id 
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Then 




du 
dr 



so the equation of motion becomes 




af 



P' 



Now let V {u) be any potential. Demanding 



-irt f 

^ " f 2 



we integrate to find 



a 



With this choice for /, the equation of motion becomes simply 

Notice that u is not necessarily a monotonic function of x so the transformation at zeros of 
V may be singular. We will not deal with such global issues here. 

In higher dimensions the regularizing transformation is complicated by the presence of 
angular momentum. Still, the general proof is similar, involving a change of both the radial 
coordinate and the time. Once again, more general potentials can be treated. To begin, we 
eliminate the angular momentum variables to reduce the problem to a single independent 
variable. The only remaining difficulty is to handle the angular momentum term in the 
radial equation. 

We end the Section with some questions: 

1. To what degree can regularizations be accomplished by canonical transformations? 
What is the relationship between regularizations and canonical transformations? 

2. What can be proved about extrema, boundedness and stability of orbits in monotonic 

central potentials bounded by various power law potentials? in monotonic central 
potentials? in arbitrary central potentials? 
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III. INEQUIVALENT LAGRANGIANS AND HAMILTONIANS 

One of the more startling influences of quantum physics on the study of classical me- 
chanics is the realization that there exist inequivalent Lagrangians determining a given set 
of classical paths. Inequivalent Lagrangians for a given problem are those whose difference 
is not a total derivative. While it is not too surprising that a given set of paths provides 
extremals for more than one functional, it is striking that some systems permit infinitely 
many Lagrangians for the same paths. There remain many open questions on this subject, 
with most of the results holding in only one dimension. 

The existence of classically inequivalent Hamiltonians is not so clear, since there are far 
more transformations preserving Hamiltonian structure than there are preserving Lagrangian 
structure. However, distinct Hamiltonians abound in quantum theory, where equivalent 
Hamiltonians may lead to distinct quantum structures ^|. If there is more than one 
Hamiltonian for a system, which one do we quantize? Furthermore, while it is clear that 
distinct Hamiltonians can lead to different quantum theories, what about the converse? Do 
there exist distinct Hamiltonian operators, H, H' with identical expectation values for all 
observables? Can distinct Hamiltonian operators have the same energy spectra? 

Here, we restrict our attention to classical questions. To begin our exploration of inequiv- 
alent Lagrangians, we describe classes of free particle Lagrangians and give some examples. 



Next we move to the theorems for 1-dim systems due to Yan, Kobussen and Leubner 
3,0,0,0) .nCuding a .„,p.e example. Then we eon.de. ,„e,.,va>e,. La,.an,ians 
in higher dimensions. Finally, we briefly examine the possibilities for inequivalent Hamilto- 
nians. 



A. General free particle Lagrangians 

There are distinct classes of Lagrangian even for free particle motion. We derive the 
classes and give an example of each, noting how Galilean invariance singles out the usual 
choice of Lagrangian. 

The most general velocity dependent free particle Lagrangian is 

S = [ f{v)dt 
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We assume the Cartesian form of the Euchdean metric, so that v — -s/Sijv^v^. The equation 
of motion is 



so the conjugate momentum 



d df 

= 

dt dv^ 



OV^ V 



is conserved. We need only solve this equation for the velocity. Separating the magnitude 
and direction, we have 

V p 

This solution is well-defined on any region in which the mapping between velocity and 
momentum is 1 — 1. This means that velocity ranges may be any of four types: v e 
(0, oo) , (0, , {vi,V2) , {vi, oo) . Which of the four types occurs depends on the singular- 
ities of fv^v. Since v^/v is a well-defined unit vector for all nonzero vi, it is /' which 
determines the range. Requiring the map from v i to pi to be single valued and finite, we 
restrict to regions of /' which are monotonic. Independent physical ranges of velocity will 
then be determined by each zero or pole of /'. In general there will be n -|- 1 such ranges 

V e {Q,Vi),{vi,V2),...,{Vn,00) 

if there are n singular points of /'. Of course it is possible that Vi = (so that on the lowest 
range, (0,^2) , zero velocity is forbidden), or f 1 = cxd so that the full range of velocities is 
allowed. Within any of these regions, the Hamiltonian formulation is well-defined and gives 
the same equations of motion as the Lagrangian formulation. 

Thus, the motion for general / may be described as follows. Picture the space of all 
velocities divided into a number of spheres centered on the origin. The radii of these spheres 
are given by the roots and poles of /'. Between any pair of spheres, momentum and velocity 
are in 1 — 1 correspondence and the motion is uniquely determined by the initial conditions. 
In these regions the velocity remains constant and the resulting motion is in a straight line. 
On spheres corresponding to zeros of /', the direction of motion is not determined by the 
equation of motion. On spheres corresponding to poles of /', no solutions exist. It is amusing 
to note that all three cases occur in practice. We now give an example of each. 
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First, consider the regular situation when /' is monotonic everywhere so the motion is 
uniquely determined to be straight lines for all possible initial velocities. The condition 
singles out the case of unconstrained Newtonian mechanics, this is the only case that is 
Galilean invariant, since Galilean boosts require the full range of velocities, v G (0, oo). 

When /' has zeros, we have situations where a complete set of initial conditions is insuf- 
ficient to determine the motion. Such a situation occurs in Lovelock, or extended, gravity, 
in which the action in d-dimensions (for d even) is of the general form 

^ = J] / R'^^ A R'^'^ A ■ ■ • A R^^ A e^' A ■ • • A e^Sabcd-efg-h 

where R"'' is the curvature 2-form, the solder form and the a^, are arbitrary constants. This 
is the most general curved spacetime gravity theory in which the field equations depend on 
no higher than second derivatives of the metric 17]. In general, the field equations depend on 
powers of the second derivatives of the metric, whereas in general relativity this dependence 
is linear. Among the solutions are certain special cases called "geometrically free" jlS^ . 
These arise as follows. For some choices of the constants ak, we may rewrite S in the form 

d/2 

■1d/20d/2 



k=0 

Suppose that for all A; = 1, . . . , n for some n in the range 2 < n < d/2, we have 

ak = a 

for some fixed value a. Then the variational equations all contain at least n — 1 factors of 

Therefore, if there is a subspace of dimension m>d — n+lof constant curvature 

for a,b = 1, . . . ,m, then the field equations are satisfied regardless of the metric on the 
complementary subspace. This is similar to the case of vanishing /', where the equation of 
motion is satisfied regardless of the direction of the velocity, 

Pi = f — = 

V 
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as long as v, but not Vi, is constant. 

Finally, suppose /' has a pole at some value Vq. Then the momentum diverges and motion 
never occurs at velocity vq- Of course, this is the case in special relativity, where the action 
of a free particle may be written as 

S = Padx'^ = — Edt + Pidx^ 



-mc^ I \l 1 dt. 



With f{v) = —mc^^yi — f^/c^, we have 



mv 



v/1 - vyc^ 

with the well known pole in momentum at f = c. 

Note that there is a complementary situation for Hamiltonians. From the Lagrangians 
for the free particle, 

S = J f{v)dt 

we have the conjugate momenta 

Pi = f P = f = 9 

V 



and Hamiltonians 



Hamilton's equations are 



H = vf'- f. 



dx^ dH ~^ ( ) 

dt dpi p 

dpi dH 
— = r = 

dt dx^ 

Once again, the constancy of the momentum is immediate. However, despite the apparent 
diversity of Hamiltonians, they are locally related by canonical transformations. The only 
distinctions are the global ones, and these exactly match those described above. 



B. Inequivalent Lagrangians 



The existence of inequivalent Lagrangians for a given physical problem seems to trace 



back to Lie 



Dirac (j^l'ISlI) was certainly well aware of the ambiguities involved in 
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passing between the 



others 
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lagrangian and Hamiltonian formulations of classical mechanics. Later, 



2q). identified certain non-canonical transformations which nonethe- 
less preserve certain Hamiltonians. A specific non-canonical transformation of the 2-dim 

n ^ 

harmonic oscillator is provided by Gelman and Saletan |2q. Bolza |22j showed that inde- 
)endent Lagrangians can give the same equations of motion and, a few years later, Kobussen 



121, Yan 



13|,[lJ|) and Okubo (|28|,[29|) independently gave systematic developments show- 
ing that an infinite number of inequivalent Lagrangians exist for 2-dim mechanical systems. 
Shortly thereafter, Leubner 16| generalized and streamlined Yan's proof to include arbitrary 
functions of two constants of motion. 

Leubner's result, the most general to date, may be stated as follows. Given any two 
constants of motion, {a, (3) , associated with the solution to a given 1-dim equation of motion, 
the solution set for any Lagrangian of the form 



L (x, X, t) 



X — V 



d{v,t) 



dv 



+ 



f {x,vo,t) — 



d{v,t) 



v=vo 



dfl 

dx H — — 
dt 



(3) 



where 



d{v,t) 



is the Jacobian, includes the same solutions locally. Notice that a and /? are 
arbitrary constants of the motion - each may be an arbitrary function of simpler constants 
such as the Hamiltonian. We argue below that in 1-dim the solution sets are locally identical, 
though 1^ provides no explicit proof. In higher dimensions there are easy counterexamples. 
We illustrate a special case of this formula, of the form 

K(x,v) 



L (x, v) 



X 



-dv 



(4) 



where K is any constant of the motion of the system. This expression is valid when the 
original Lagrangian has no explicit time dependence. Following Okubo |29^, we prove that 
eq.(jl} leads to the constancy of K. The result follows immediately from the Euler-Lagrange 
expression for L : 



K (x,v) , K (x,x] 
— dv + X- 



v" 



d dL dL d 

dt dx dx dt 

X dK (x, x) dK (x, x) 

X dx dx 
1 dK (x, x) 

X 



x^ 



— X 



1 dK (x, v) 
f 2 dx 



dv 



dt 
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Therefore, the Euler-Lagrange equation holds if and only if K ( constant of the 

motion. 

The uniqueness in 1-dim follows from the fact that a single constant of the motion is 
sufficient to determine the solution curves up to the initial point. The uniqueness also 
depends on there being only a single Euler-Lagrange equation. These observations lead us 
to a higher dimensional result below. 

It is interesting to notice that we can derive this form for L, but with K replaced by the 
Hamiltonian, by inverting the usual expression, 

ox 

for the Hamiltonian in terms of the Lagrangian. First, rewrite the right side as: 

OX ox \x 

Now, dividing by x and integrating (regarding if as a function of the velocity) we find: 

L = X f — ^"^^ '^^ dv 
J 

The remarkable fact is that the Hamiltonian may be replaced by any constant of the motion 
in this expression. Conversely, suppose we begin with the Lagrangian in terms of an arbitrary 
constant of motion, K, according to eq.(jll), 

r K (x,v) 

L(x,v)=x / — dv 

J 

Then constructing the conserved Hamiltonian, 



OX \ J J J V"^ 



V^ X J I v^ 



= K {x, x) 

we arrive at the chosen constant of motion! This proves the Gelman-Saletan-Currie con- 
motion gives rise to a possible 



jecture j26[: any nontrivial time-independent constant o 

Hamiltonian. Proofs of the conjecture are due to Yan and Leubner 3| 
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The conjugate momentum to L constructed according to eq.(jll) is 

_ dL 



^ dx 



dx \ I f 2 



^ K (x,v) , K (x,x) 
dv H : 



f 2 X 



Of course, ii K = ^mx^ + V, both H and p reduce to the usual expressions. 

The simple harmonic oscillator is sufficient to illustrate the method (j3|,j3^). Since the 
Hamiltonian, H = ^mv'^ + is a constant of the motion so is H"^, so we write 



1 r 1 



4: V 



,2 



(m^f ^ + 2kmv'^x^ + fc^x*^) dv 



1 1 1 

^ 2 , 1 -22 ^12 4 

= — m X H — kmx x — -k x . 
12 2 4 

The Euler-Lagrange equation resulting from L is 





d dL 


dL 


dt dec 


dx 


d /I 


2 • 

m X ■ 


di \3 



— {kmx^x — k^x^) 



= (mx + kx) (mx^ + A;x^) . 

Either of the two factors may be zero. Setting the ffist to zero is gives the usual equation 
for the oscillator, while setting the second to zero we find the same solutions in exponential 
form: 

The conjugate momentum and Hamiltonian are: 

H {x,p) = {x,p) = - (jn^x^ + 2kmx^x^ + fc^x'^) 

~ 2 • S 7 -2 

p = — — = -m X + kmxx . 
dx 3 

While it is possible to solve the cubic equation to find x {p) , and then substitute to find 
H (x, p) as an explicit function of p, it is clear that the resulting expression is not of the 
same form as the original Hamiltonian. There remains the question of whether this effect 
could be achieved by a time-independent canonical transformation. The transformation of 
the momentum, 

p = h kpx'^ 

6m 
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is part of a canonical transformation, given by 



p — h kpx"^ 

3m 

However, this does not simplify the form of the Hamiltonian. We find: 




H {x,p) = ^ (^^'^''^ + m^p^e^'^* + ^mpe-^^^^ 



and the resulting Hamiltonian equations of motion are not transparent. 

There does exist, of course, a time-dependent canonical transformation relating the two 
Hamiltonians. The systems are nonetheless distinct globally, since the cubic relationship 
between momentum and velocity limits the allowed ranges of the variables for the higher 
order Hamiltonian. It would be interesting to know if there is a unique Hamiltonian for 
which p {v) is 1 — 1. 



Inequivalent Lagrangians have been defined as Lagrangians which lead to the same equa- 
tions of motion but differ by more than a total derivative. For the simple case above, the 
cubic order equation of motion factors into the energy times the usual equation of motion, 
and setting either factor to zero gives the usual solution and only the usual solution. How- 
ever, is this true in general? The Yan-Leubner proof shows that the new Lagrangian has 
the same solutions, but how do wc know that none of the higher order Lagrangians intro- 
duces spurious solutions? The proofs do not address this question explicitly. If some of 
these Lagrangians introduce extra solutions, then they are not really describing the same 
motions. 

Suppose, for some time-independent Hamiltonian we write 



where a is any constant of the motion. Then we know that the Euler-Lagrange equation is 
satisfied by the usual equation of motion. But what is the Euler-Lagrange equation? We 
have shown that 



1. Are inequivalent Lagrangians equivalent? 




d dL dL ldK{x,x) 1 ,da{x,x) 



dt dx dx X dt X dt 
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Setting this to zero, we have two types of solution 

f («) = 
da Q 
'dt ~ 

If spurious solutions could arise from motions with /' = 0, those motions would have to 
stay at the critical point, say, of /. But this means that a = aQ remains constant. 
Therefore, the only way to introduce spurious solutions is ii da/dt = has solutions beyond 
the usual solutions. This may not be possible in one dimension. Finally, the inverse of the 
equation a {x, t) = ao may not exist at critical points, so the theorem must refer only to 
local equivalence of the solutions for inequivalent Lagrangians. 



C. Inequivalent Lagrangians in higher dimensions 

It is of interest to extend the results on inequivalent systems to higher dimension. Pre- 
sumably, the theorems generalize in some way, but while one dimensional problems may be 
preferable "for simplicity" [16|, this restricted case has many special properties that may 
not generalize. In any case, the method of proof of the Kobussen-Yan-Leubner theorem does 
not immediately generalize. 

For 1-dim classical mechanics, there are only two independent constants of motion. The 
Kobussen-Yan-Leubner theorem, eq.Q, makes use of one or both to characterize the La- 
grangian and, as noted above, one constant can completely determine the paths motion in 
1-dim. The remaining constant is required only to specify the initial point of the motion. 
This leads to a simple conjecture for higher dimensions, namely, that the paths are in general 
determined by n of the 2n constants of motion. This is because n of the constants specify 
the initial position, while the remaining constants determine the paths. 

We make these comments concrete with two examples. First, consider again the free 
particle in n-dim. The usual Hamiltonian is 




and we immediately find that a complete solution is characterized by the initial components 
of the momentum, poi and the initial position, xqi. Clearly, knowledge of the momenta is 
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necessary and sufficient to determine a set of flows. If we consider inequivalent Lagrangians 



where 



then the momenta 



comprise a set of first integrals of the motion. Inverting for the velocity 

v' = v' (pio) 

fixes the fiow without fixing the initial point. 

In general we will need at least this same set of relations, — (pio) , to determine the 
fiow, though the generic case will involve n relations depending on 2n constants: 

v' ^v' {pio,xi) . 

Notice that fewer relations does not determine the fiow even for free motion in two dimen- 
sions. Thus, knowing only 

POx 
Vx^ — 

m 

leaves the motion in the y direction fully arbitrary. 

In an arbitrary number of dimensions, we find that expression for the energy in terms 
of the Lagrangian is still integrable as in the 1-dim case above, as long v — ^A'^. If the 
Lagrangian does not depend explicitly on time, then energy is conserved. Then, letting 
6'^ — xY^, we can still write the Lagrangian as an integral over Hamiltonian: 

L(^,v,e,j=vJ ^ d^ + f[^,0,j 

where / ^x, is now necessary in order for L to satisfy the Euler-Lagrange equations. 
The integral term of this expression satisfies one of the Euler-Lagrangc equations. If we 
now define a new Lagrangian by replacing H by an arbitrary, time-independent constant of 
the motion, a ^x, v, , 

L^v j ^ ^ dc + f{x,e) 
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then the new Lagrangian, L, still satisfies the same Euler-Lagrange equation 

J d dL dL\ 
\dtdx' dx'l 

We conjecture that for a suitable choice of /, L provides an inequivalent Lagrangian, thereby 
providing one of the n relations required to specify the fiow. 



D. Inequivalent Hamiltonians 

The question of inequivalent Hamiltonians is quite distinct from that of inequivalent La- 
grangians, because the group of allowed transformations is much larger. Indeed, Hamilton- 
Jacobi theory shows that any Hamiltonian may be made trivial by a canonical transforma- 
tion. This also means that any two Hamiltonians are related locally by a time-dependent 
canonical transformation. At least in this sense, all Hamiltonians with the same number of 
degrees of freeedom are locally equivalent. 

In light of this, the first question to be answered is the following. Since the properties 
of canonical transformations are defined by the demand that they change the action by 
no more than a total derivative (for example, in the derivation of the properties of gener- 
ating functions), how can two Hamiltonians be locally equivalent while the corresponding 
Lagrangians are inequivalent? The answer is subtle. Canonical transformations are defined 
in such a way as to leave the Hilbert one-form 

dx^ 

L (pi, x\ q\ t) dt = Pi-^dt - H {pj, q^) dt 
changed by no more than an exact form. But this L is not quite the Lagrangian, since 

L = L {pi,x\q\t) 

while a true Lagrangian is a function of and only. Since the Lagrangian formalism is 
invariant only under coordinate diffeomorphisms, = a;* {q-^,t) , canonical transformations 
involving both pi and are not expected to preserve it. 

Despite Hamilton- Jacobi theory, there are ways to define a notion of inequivalent Hamil- 
tonians. First, if we restrict to time-independent canonical transformations, there exist 
distinct Hamiltonian systems related by diffeomorphisms, H{q,Tr) — H[x {q,Tr) ,p {q,Tr)]. 
This allows substantial variation in the functional form of the Hamiltonian and it may be 
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difficult to determine whether two Hamiltonians are related in this way. Second, we may 

define canonoid transformations, defined as preserving the canonical structure for one or 
I I 

more Hamilton [26]. Presumably, these are related to symmetries of particular systems. Or, 
third, we may quantize the system and ask whether the quantum systems are equivalent. 

These considerations point to a heirarchy of classifications of Hamiltonian systems, equiv- 
alent up to some set of transformations. It would be useful to know the exact set of trans- 
formations under which a given set of phase space paths is invariant. We know that the set 
is smaller than time-dependent transformations and larger than canonical transformations, 
and seems likely that the answer depends on the class of curves in some way. Since it is 
really the solution curves that define equivalence, it is clear that some systems will have 
more symmetry than others. 

Even if we had a clear characterization of Hamiltonians for a given system, it is not clear 
that we know what the system is. For example, suppose a given system admits a time- 
independent canonical tranformation taking the Hamiltonian to a constant. Such a system 
exists if we can find nontrivial solutions to the time-independent Hamilton- Jacobi equation 

— VS-VS + V = 

and this is surely possible for some systems. But this means that one of the equivalent 
formulations of the problem describes straight-line motion. Clearly, finding the solution 
curves is not enough to describe such a system - we must also keep track of the sequence 
of transformations we used to trivialize that solution. Thus, specifying a classical system 
requires some statement about the correspondence between phase space coordinates and 
measurements in some physical system. Of course, there are many equivalent ways to set up 
such a correspondence, but at least one must be specified and the subsequent transformations 
tracked. 



IV. CONSTANTS OF MOTION 



Recent decades have seen interesting new techniques and revivals of known results for 
symmetries ( 3^,j3^). Some of these have to do with the Kepler problem. The best-known 
rediscovery concerning the Kepler problem is that in addition to the energy, E and angular 
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momentum, 



^ 1.2 a 

E = -mx 

2 r 

L = r X p (5) 

the Laplace- Runge-Lenz vector (0|, 0|, 0], 0]) is conserved. We define the Laplace- 
Runge-Lenz vector by 

A = p X L — mar. 

Geometrically, A points in the direction of the periapsis, and may therefore be thought of 
as giving specifying the orientation of the orbit within the orbital plane. Keplerian orbits 
can be described completely in terms of six initial conditions, and since one of these is the 
initial position on a given ellipse, only five remain among the energy, angular momentum 
and Laplace- Runge-Lenz vector j38|. Two constraints - the orthogonality of A and L, and a 
relationship between the magnitudes A, L and E - give the correct count. Of course, these 
three quantities are not the only set of constants we can choose. A number of fairly recent 
3^;|i^;|il|) have identified a simpler conserved vector quantity, which (lacking 



authors 



evidence for an earlier reference) we will call the Hamilton vector 4^. It is given by 



u = v--c, 

and may be used together with the energy and angular momentum as a complete set of 
constants. Apparently this vector was well-known in the 19th century, then dropped from 
texts j^. Its time constancy is a direct consequence of the force law since, for an arbitrary 
central force /(r), 

duj _ dvj ad0i _ f (r) Xj ^ a dip Xj 
dt dt L dt mr mr'^Lp dt r 



a 



X. 



-fir) — (7) 
J mr 

and dui/dt = precisely when /(r) is given by an inverse square law . Notice that it is the 
balance between the radial dependences of angular momentum and the force that allow this 
characterization. We give a derivation of the Hamilton vector below. 

The Laplace-Runge-Lenz vector, the Hamilton vector and the angular momentum are 
related by 

A = mu X L 
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and 



L X A = mL X (u X L) 



(8) 
(9) 



where we have used the fact that u hes in the plane of the orbit and is therefore perpendicular 
to L. 

It might be of interest to study more general central force problems using time- dependent 
versions of the Laplace-Runge-Lenz and Hamiltonian vectors. While no longer conserved, 
they are more geometrical than the usual polar coordinates. This could result in some 
simplification. It would also be of interest to know if these vectors correspond to symmetries 
- perhaps they reflect symmetries of the corresponding phase space, or, since the bound state 
Kepler problem may be embedded with an 50(4) symmetry, perhaps they are part of that 
symmetry. It seems likely that the Kepler problem has an even larger symmetry - perhaps 
5*0(4, 1) - since the open orbits have 50(3, 1) symmetry. This might be explored by writing 
the Kepler problem in terms of Spin (4, 1) conformal spinors. 

Munoz j2| shows that the Hamiltion vector leads to an easy derivation of the equation of 
motion. Indeed, let the perihelion of the orbit occur at time t = on the x-axis so that the 
velocity is given by 

V = Vo(p 

Then u = uj, where the unit vector in the y-direction gives the initial direction of 0. Dotting 
u with if we have 

u ■ (fi = v ■ (fi — a/L 

u COS0 = rip — a/ L 



or replacing = L/{mr'^ 



1 mu am 
- = —oosip + — 

L?' /ma 



r 
r 



(10) 



1 + (Luo/a) cos(p 
as usual. 

Tjiang and Sutanto fl| describe a straightforward way to identify constants of the motion 
arising from the vanishing of 



df . , df 
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based on a well-known theorem on the solution of differential equations. The theorem states 
that any equation of the form 

i * 

has the general solution given by 

/ = <l>(Mi,...,Mfc) 

where k <n and the (xi, . . . , Xn, /) are solutions to 

dxi dx2 dxn df 

Pi P2 Pn R 

Applied to eq. ljllll . the theorem implies that the functions Ui are all constants of the motion. 
Moreover, we can compute the possible constants of motion by solving the equations 

dxi dXn dpi dXn 



dH \ I dH \ (_dH_] {—ML 

dpi I \ dpn J \ dxi J \ dxn 



dt 



To illustrate the method and at the same time derive the Hamilton vector as a constant 
of Keplerian motion, we apply the technique to the Kepler problem. First, the Hamiltonian 
is given, in any dimension n > 2, by 

2m r 

so we must solve the equations 

— dxi = ■ ■ ■ = — dXn = dpi = ■ ■ ■ = dpn = dt (121 

Pi Pn OtXi aXn 

First, for each consider the equations of the form 

m , 



— dxi = dpi 

Pi (XXi 

aXi 1 

-dXi H Pidpi = 0. 



^3 



Summing over i we have the constancy of the Hamiltonian: 

dH = d(T^-^)=0 
Next, consider the equations among the dxi. For any pair of these {i 7^ j)we have 
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so that 

= pjdxi — Pidxj (13) 
= d {pjXi - piXj) - {dpjXi - dpiXj) (14) 

We may replace the momentum differentials using the Tjiang-Sutanto equations, eqs. lfT^ to 
write 



dpjXi — dpiXj = dxi H -dxi = 



axj mxi , axi mx 
dxi H 5 

so we have conservation of all of the components 



iWjj XiPj XjPi 

of the angular momentum. Of course, in 3-dimensions we may use the Levi-Civita tensor to 
write this as 

Li = ^Sijk {xiPj - XjPi) = [x X p] . 

In Sec. 2 we showed from the constancy of Mij that the motion remains in a fixed plane for 
all time. 

Finally, we study the additional constants of motion arising from equations of the form 

m 

— dxi = dpi 

Pi aXi 

Since Pi = mdxi/dt, this may be written as 

TflT 

dt + dvi = (15) 

^dt + dvi = (16) 



-2, 



Now, using L = mr Lp where L is the magnitude of the angular momentum, L = 
a/X] Mij Mij, and (p is in the plane of the orbit, we have 

dvi + ^—dif = (17) 
L r 

+ ^ sin = (18) 

and similarly 

d (^vi — ^ cos (fj = 
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Adding these with unit vectors and noting that v = vii + V2J comprises the entire velocity 
vector, estabhshes the constancy of the Hamilton vector: 

u = v--c, 

Thus, {if, Mjj, Ui} is a complete set of constants of the motion for the n-dim Kepler problem. 
The solution, eq. |TO|) . follows immediately. 



V. HIGHER ORDER EQUATIONS OF MOTION 

A. Euler-Lagrange and Hamiltonian systems of arbitrary order 

There generalization of the Euler-Lagrange equation to systems for which the Lagrangian 
depends on higher than second derivatives of the postion is immediate and well-known. If 

J-J J-J ( X ^ Ob ^ Ob ^ • • • ^ OC , t) the resulting variation leads to 

jk=0 

This generalized Euler-Lagrange equation is generically of order 2n. Such systems are, of 
course, allowed within even the Newtonian formulation. A simple electronic circuit and 
feedback mechanism can easily drive a motor in a way dependent upon rates of change of 
the acceleration. Moreover, there are systems of genuine physical and mathematical interest 
which require higher order differential equations for their description. After deriving a few 
general results for higher order systems, we look at one example in detail - Bohmian quantum 
mechanics. 

Returning to generalised Euler-Lagrange systems, suppose L is independent of time. Then 
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But 



so 



dxi^) dt\ dxi^)) dtdxi^) ^ ' 

= ± (x^^^^\ - - (x^^-^^-^\ r2oi 

dt\ dxi^)) dt\ dtdxi^)) ^ ' 

^ dt^dxW ^ ^ 

; (22) 

- \^ f^-]"^ A ( Ak-m)^_^J^ _ (1) (_\k^l (no\ 

' dt\ dt^^dxW) ^ ' dt^^dxO^) ^ ^ 



dL _ d s^sr^ f (k^rn) d"" \ ^ , d'' dL 

dt dt^^^ ' \ dt^ dxi^) ^ 



, dt^dx^^) 

k=0 m=0 ^ ' k=0 



Using the equation of motion, 

E, . ^ d^ OL 

k=Q 

the final sum vanishes and we have the conserved energy 

n k—l 



I — n ^ — n \ / 



fc=0 m=0 

The n = 3 case of this result is given in j43| and elsewhere. 

Directly from the generalized equation we see immediately that if a coordinate x is cyclic, 
(i.e. dL/dx = 0), we still get a conserved momentum, 

E/ \m+l d"^ 9L 



dV^ (9x("+^) 

m=0 



This follows from 



dt^ dxC") 

k=0 



k=l 

m=0 
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With higher order Lagrangians, there are additional possibihties. Suppose the lowest m < n 
partials of L vanish: 

dL 

0, fc = 0, l,...,m-l 



Then the sum in the field equation starts at m, and extracting m time derivatives 



k=m 



SO that the momentum 

= ^ I V (-) 



k=m 



is conserved. Integrating m — 1 more times, 



m—l 



E W = E (-)' 



k=0 k=m 

where we now have m constants, p^. 

Higher order systems also permit a Hamiltonian formulation. Let n — 2m — 1 be any 
odd integer. We divide the time derivatives of x into even and odd order, and replace the 
odd time derivatives with conjugate momenta. For even n, Hamilton's equations will be 
supplemented by one additional Euler- Lagrange equation. Thus, let 

dL 

Pk = 



for A; = 0, 1, . . . , m — 1. The Legendre transformation is employed in the usual way to express 
the Hamiltonian in terms of and the Pk, 

dyk 



H (yk,Pk) = ^Pk-^ - L 



k=l 



As usual, H is independent of the odd accelerations, x^'^^~^^\ since 

dH dH dL 

Pk - (nl.,U = 



The variation of the Lagrangian with respect to yk and pk is straightforward: 

= 6S 



Y^PkSy + J2 y^Pk - E d:;rjyk - E d^jpk 

k=l k=l k=l ^k i,=l ^k 
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Integrating by parts we find 2m = n + 1 first order equations: 

dyk _ dH 

dt dpk 
dpk dH 

dt dyk 

These are the generahzed Hamilton's equations. 

Naturally, higher order equations require more initial data than we usually have to specify 
to determine the motion of a classical system, so their occurrence is somewhat rare. But 
ultimately, any restriction to second order equations in classical physics is phenomenological, 
depending principally on the success of second order models for fitting measurements. If 
we take quantum physics into account, higher order equations in field theory may introduce 
ghosts or other undesirable features. 

Nonetheless, there are situations where higher order equations are justified. We briefly 
discuss one of these below, Bohmian quantum mechanics. The KdV equation provides an 
example of an integrable system, having an infinity of independent constants of motion. It 
is presented in an Appendix, as is the approach to chaos. 



B. Bohmian quantum mechanics 

A central theme of the Bohmian approach to quantum mechanics is to give it a form 
which may be interpreted classically The first step is to replace the complex 

wave function by pair of real valued functions. This is accomplished as follows. Let 

where A and 5* are real valued functions. Substituting into the Schrodinger equation, 

and separating the real and imaginary parts give two equations: 

^VS .VS + V + ^ = i^^V'A (28) 
2m dt 2mA ^ ^ 

dA 1 ^ ^ 1 

— + — V5 • VA + — AV^S = (29) 
dt m 2m 

We refer to this system as the Bohm equations. The first equation is the Hamilton- Jacobi 
equation with an additional term. This equation is frequently used to show how the classical 
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limit emerges from quantum mechanics when h ^ 0, but here we want to exactly replicate 
the content of the quantum theory while maintaining a classical viewpoint. Multiplying the 
second equation, eq. pi?|) . by A, it may be rewritten as 

^ + -V ■ (a'Vs) = (30) 

Defining 

P = A2, v = -V^ 
m 

the equation becomes the continuity equation for a current J = pv, 

This is the usual conserved probability current of quantum mechanics, cast in classical guise. 
Alternatively, we may view eq.(j2Hl) as a wave equation for A : 



h',, , f 1 , , ds\ 



+ — VS ■V^ + V^+— M = 

2m \2m ot J 

This is just a diffusion equation with a horrible potential. 

In the one dimensional, stationary case we can reduce the Bohm equations to a single, 
higher order differential equation. In 1-dim, with the stationary conditions 

as ^ dA 

the continuity equation, eq. ljHUj) is simply 

- (A'S'Y = 
m 

which integrates immediately to give 

^ - S' 

Now substituting this result into eq.(j2HI) together with dS/dt = —E, results in a third-order, 
nonlinear equation 



2m 



S'" 3/3"^^ 



2S' 4 V S' 



+ 7^ (Sf + V-E = (31) 
2m 



The interesting point is that this equation is rigorously equivalent to the 1-dim stationary 
state Schrodinger equation. The downside is that we have handled only the 1-dim, stationary 
case. 
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This equation lacks the simplicity of the Schrodinger equation. For example, suppose we 
solve 



2m 



S' 



2S' 



— is'f + y - £; = 



to find a solution S {x,E) . Unfortunately, the nonlinearity means that we cannot take a 
superposition of stationary states to get a general time-dependent solution. 

Consider the time-dependent case further. We can at least find a time dependent solution 
to the (linear!) second equation: 



^ + -S'A' + —AS" = 
at m 2m 



Direct integration shows that 



m 



dx 
'S' 



is a solution for any function G. The amplitude therefore propagates with fixed spatial form. 
It would be of interest to know if this type of solution for A generalizes to higher dimensions. 
Now we may substitute this into the first equation, 

1 



^ + — {sf + v ^ 

^ 2m 2m^^^^_^j|^ 




^Zcit-mlt 



Perhaps an appropriate separation can solve this equation as well. 



1. Bohmian Lagrangian 

The 1-dim, stationary state Bohm equation also follows from the variation of the usual 
Schrodinger action, reduced by the solution for A. Starting from the Schrodinger action, 

we substitute the polar expression t/j — -^ei''^^^^"^*! for the wave function, the action 
becomes 

r {S'f 1 v-E 

J 2m 4 {S')^ 2m S' 
Varying, the equation of motion is found to be 

2S"' h'' 3S"S" ff 3 (S'f 1 V-E 
2m4:{S'f ^ 2m2{S'f ~ '^rn 4 {S'f ^2m~ (S'f 



d 

dx 
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The term in brackets is a constant. Choosing the constant to be 1/m correctly reproduces 
the stationary Bohm equation: 

AmS' 8m\S'J 2m ^ ' 
We could have avoided the need to pick this integration constant by taking the action to be 



This form differs from the previous one only by the integral of a total derivative. 
Writing the equation of motion in the form, 

= -^S'S'" + — {S"f - — {S'f -{V-E) {S'f (32) 
4m 8m ^ ^ 2m ^ ^ ^ ' ^ ' 

it has been observed (j43],|3|) that there are solutions with S' = 0. It is not difficult to see 
that these should not be considered to be the physical solutions. First, we arrived at eq.(jH21) 
by the substitution 



which is singular for constant S. Since wave functions with divergent amplitude are not 
considered physical, such points require a closer examination at the very least. Second, we 
can see that these are clearly not the physical solutions because they are independent of the 
potential V. 

To eliminate the spurious solutions, while indicating a method of solution, we make a 
simple change of variable. Notice that the nonlinear terms in the Bohm equation, eq. ()3ip . 
may be rewritten as 




V-E 
S' 





Substituting, the equation of motion becomes 




Now let r 





2mr^ 



1 



+ {V -E)r = 



Rearranging, 



m?/h^ 2m? 



iy -E)r 



mr — 



mr^ h? 
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the form suggest that we think of the independent variable x as time. Replacing x — t, we 
interpret r {t) as a radial coordinate and \^ as a time-dependent potential V (t) . Then we 
have 

mr = -—[V{t)-E]r 

and we recognize the isotropic 2-dim harmonic oscillator with angular momentum m/h with 
a time-dependent spring strength, 

k{t) = -^[V{t)-E] 

This clearly has bound state solutions for suitable energies and potentials. Notice that 
solutions with 5" = correspond to infinite radial coordinate, so bound state solutions 
automatically avoid this spurious case. 

It is suggestive that the introduction of a time-dependent spring constant into the 
isotropic oscillator can lead to parametric resonance It would be of interest to find the 
relationship of such parametric resonances to the eigenmodes of the corresponding quantum 
problem. 



VI. GAUGING NEWTON'S LAW 



One surprising new result in classical mechanics is that both the Lagrangian and Hamil- 
tonian formulations of Newton's laws may be derived as gauge theories of Newton's second 
law To see how this comes about, and to understand the symmetries involved, 

we digress a moment to consider the essential elements of a physical theory. In particular, 
we want to distinguish two features: dynamical laws and measurement theory. 

The distinction between these is easy to see. For example, in quantum mechanics the 
dynamical law is the Schrodinger equation 

Hip - ih— 

ot 

which governs the time evolution of the wave function, ip. The measurement theory is what 
establishes the correspondence between calculations and measurable numbers. One of the 
chief elements of quantum measurement theory is therefore the Hermitian inner product on 
Hilbert space: 

Jv 
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As a second example, consider Newtonian mechanics. The dynamical variable for a particle 
is the position vector, x, and its motion is governed by the second law: 



while the inner product allows us to extract measurable magnitudes 

(u, v) = U • V 

We are intested in the differences in the symmetries of dynamical laws and measurement. 
Generally, the differential or other equation governing dynamical evolution is invariant under 
some global symmetry. In contrast to this, a metric, an invariant product or some other real- 
valued mapping to measurable quantities is often invariant under a group of diffeomorphisms. 
Whatever the symmetries, it often occurs that the symmetry of dynamical evolution and the 
symmetry of the measurement theory are different. Which is the symmetry of the system? 

For Newtonian measurements, the inner products allow local transformations and there- 
fore have the larger symmetry. It makes sense to try to extend the symmetry of the dy- 
namical law to agree with that of the measurement theory. Fortunately, there are standard 
techniques for accomplishing this extension - the methods of gauge theory. 

Gauging takes a global symmetry, that is, a symmetry that is independent of position and 
time, and extends it to a local symmetry, i.e., one that may be different at different positions. 
We systematically extend to a local symmetry by introducing a connection: a one-form field 
vahicd in the Lie algebra of the symmetry wc wish to gauge. Added to the usual partial 
derivative, the connection subtracts back out the extra terms arising from differentiating 
the local symmetry. The most familiar example is general relativity, in which the Christoffel 
connection, F^^ added to partial derivative makes the derivative covariant with respect to 
general coordinate transformations. The U{1) gauge theory of electromagnetism is also 
familiar. In this case, the vector potential provides the connection. 

What arc the symmetries of the Newtonian dynamical and measurement theories? There 
is more than one answer. The dynamical law is invariant under the Galilean group, G, 
consisting of rotations, translations, Galilean boosts and time translations. It is possible 
to extend the rotations to general linear transformations and still leave the second law 
invariant. For the measurement theory, the Euclidean line element is invariant under the set 
of 3-dimensional rotations and translations, I SO (3). This is called the Euclidean group. If 



35 



we regard Euclidean 3-space as a manifold instead of a vector space, these transformations 
may be local. Furthermore, recognizing that we only actually measure dimensionless ratios 
(for example, the ratio of the height of a tree to the length of a meter stick), we can require 
invariance of ratios of line elements. This gives the conformal group, 50(4, 1). 

We will consider two of the possible gaugings of Newton's second law - Euclidean and 
conformal. 



A. Euclidean gauge theory of Newton's second law 

The Euclidean gauging of Newton's second law leads to Lagrangian mechanics. This 
is not a particularly surprising result. However, while it is well-known that Lagrangian 
mechanics provides a form of Newton's second law valid for "generalized coordinates," the 
construction by gauging has advantages: (1) it arrives at the Lagrangian formulation in a 
systematic way, (2) gauging displays explicitly the meaning of generalized coordinates, and 
(3) it illustrates the general techniques used for the gauging of conformal symmetry below. 
As we shall see, the results of conformal gauging in the next subsection are unexpected. 

Proceeding, first recall that the transformations of the Euclidean group IS0{3) include 
three rotations and three translations, and the Lie algeba has corresponding generators. 
Gauging therefore gives us two sets of 1-form gauge fields: 

1. Three translational gauge fields, comprising the dreibein, e*. 

2. Three rotational gauge fields, the S0{3) spin connection, Uj, antisymmetric under the 
interchange of indices. 

The gaugin g pr oceeds just as when we gauge the Poincare group to develop Riemannian 



geometry (pi|, In Poincare gauging, the vierbein, e", is identified with an orthonormal 
frame field on a 4-dim Riemannian manifold and the spin connection permits the use of 
local Lorentz transformations. For the present Euclidean case, the dreibein, e*, is identified 
with an orthonormal basis of a 3-dim manifold and the 5*0(3) spin connection, cuj, gives local 
rotational symmetry. The pair (e*,Ci;*) is equivalent to the metric and general coordinate 
connection, {gmn, P™) . 
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The connection forms must satisfy the Lie algebra relations of the symmetry group, as 
encoded in the Maurer-Cartan structure equations: 

The solution of these is simple since we do not include curvature. The first equation is solved 
by the pure gauge form of the connection, 

where A*- (x) is a local rotation matrix. This means that there exists a choice of frames (say, 
Aj — constant) in which the spin connection is zero. Choosing this frame, the equation for 
the dreibein is satisfied by setting 

= dx' 

Prom this we see that the equations describe Euclidean 3-space. Using the spin connection 
we can define a derivative operator which is covariant with respect to local rotations. If 
we cast the same equations in terms of a coordinate basis using the metric and Christoffel 

connection, {gmn, Tr") ' the derivative is covariant with respect to general coordinate changes, 
or diffcomorphisms. 

We may find the new dynamical law using a variational principle. Using the coordinate 
metric, 

9mn ^ ^ij 

we choose the squared norm of the velocity vector, plus a function of the coordinates to 
provide a source for the motion: 

S = j[gmnv'^v^ + (l>{xn]dt 

Because we have local symmetry, we can write the same thing in any coordinates. Notice 
that there is always some arbitrariness in the gauging procedure at this point. There are 
two properties we demand of this variational principle. First, it must be invariant under 
the local symmetry group. Second, we require the restriction of the new dynamical law 
to the original symmetry to reproduce the original law. The action S above satisfies these 
requirements. 
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Varying S, we find the new form of Newton's law, 

Dv"" _ dcf) 

where the covariant derivative of transforms as a vector under local rotations. When 
= 0, this is the geodesic equation. Since the space is Euclidean, the geodesies are straight 
lines. The class of straight fines 

is equivalent to the class of Newtonian inertial reference frames. 

Writing V — —^ip lor & potential V, we see that forces produce deviations from geodesic 
motion. This is the Lagrangian formulation of mechanics. Note that we get the same 
equation of motion if we substitute the Lagrangian in the form 

into the usual Euler-Lagrange equation. The general coordinate invariance ( "use of general- 
ized coordinates") is, of course, one of the main reasons for the use of Lagrangian methods. 
The present approach, while principally intended to pave the way for the conformal gauging 
below, does have the advantage of systematically showing that the class of generalized coor- 
dinates is just the diffeomorphism group. In the usual formulation, this conclusion follows 
from the coordinate invariance of the action. 



B. Conformal gauge theory of Newton's second law 

We now repeat the gauging process, but this time use the full conformal symmetry. The 
conformal group (for compactified 3-dim Euclidean space) contains ten transformations: 

1. 3 rotations 

2. 3 translations 

3. 1 dilatation 

4. 3 special conformal transformations 



38 



The first two sets of transformations reproduce the Euchdean group. Dilatations just 
rescale all lengths by a factor, while special conformal transformations are translations in 
inverse coordinates. 

These global transformations preserve the Euclidean line element up to an overall multi- 
ple. As it stands, Newton's second law is not invariant under even the global form of these 
transformations - the special conformal transformations do not leave the law unchanged 
because they do not act linearly on Euclidean 3-space. This is easy to fix: we introduce 
a very limited covariant derivative with a connection specific to global special conformal 
transformations. It is unusual to require a connection in a dynamical law before gauging, 
but nothing forbids it and it gives us an equation with the symmetry we wish to gauge. 

Now consider Newton's law, modified just enough to let us perform all 10 global conformal 
transformations. Make those 10 global transformations local. There is more than one way 
to do this, but so far only one appears to be interesting - the biconformal gauging described 
below. 

There will now be ten gauge fields: 

1. The dreibein, e* 

2. The (antisymmetric) 5*0(3) spin connection, a;] 

3. The Weyl vector, W. 

4. The co-dreibein, fj, from special conformal transformations 

The local rotations, gauged by the spin connection, are as expected and we add local 
dilatations gauged by the Weyl vector. These allow general coordinate invariance and scale 
invariance. Employing the biconformal technique, we interpret (e*,fj) as an orthonormal 
frame field of a six dimensional space. 

These gauge fields must satisfy the Maurer-Cartan structure equations: 

duj] = uo^uol + e'ij - e^r 

de^ = ^uj] + We^ 

df, = a;ff,- + f,W 
dW = 2e% 



39 



This is just the conformal Lie algebra in a dual basis. Once again the equations are easily 
solved. The solution reveals a symplectic form, 

e = e% 

The six dimensional space therefore has a similar structure to a one particle phase space. The 
units of the coordinates of this 6-dim space are not all the same. Three are correct for posti- 
tion (x*, length) while the remaining three are geometric units for momentum 1/ length). 
Note that the conversion of units of momentum to units of inverse length may be accom- 
plished using any conventional dimensional standards, e.g., meters, seconds and kilograms. 
It also follows from the solution that the Weyl vector is given by 

W = -yM 

To find the new dynamical law we again write an action. Since the geometry is like phase 
space, the paths will not be anything like geodesies, so path length will not work. Instead, 
we have a new feature - the Weyl vector - that comes from the dilatations. We will base 
our dynamical law on the geometric interpretation of this vector field. We digress briefly to 
explore its properties. 

It follows from the nature of conformal geometry that the integral of the Weyl vector 
along any path gives the relative physical size change along that path: 

I = lo exp J Wiv'^dt = lo exp J W 

This means that magnitudes are not preserved - initially identical rods transported along 
different curves might be different sizes when they are returned together and compared. 
This possibility is the price we pay for the freedom to make local scale transformations, just 
as in a Riemannian geometry vectors may rotate even under "parallel" transport. We will 
return to this point below. 

We take the action to be the integral of the Weyl vector. Then the physical paths will be 
paths of extremal size change. Notice that, while the exponential above is gauge dependent, 
its variation is not. Indeed, it is worth noting that the gauge freedom of the Weyl vector 
agrees exactly with the freedom to add a total derivative to a Lagrangian. 
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Once again we add a function to provide a source, 

S= f (W-Y + (h)dt 



It is interesting that such a function is provided automatically in the relativistic version of 
this gauging, as the time component of the Weyl vector. 
Now vary the action. There are six first-order equations: 

(ix* d(f) 

dt dyi 
dtji _ d(j) 

dt dx^ 

If we identify (p with the Hamiltonian, these are Hamilton's equations. Therefore, the 
gauge theory of Newton's second law with respect to the conformal grouip is Hamiltonian 
mechanics. 

There are a couple of points to be clarified. First, the multiparticle case works even though 
a single Weyl vector must account for the Hamiltonian and momentum of each particle as 
long as we assume that two particles never occupy exactly the same space. This is consistent 
with the usual requirements of Newtonian mechanics, by which matter is impenetrable. 

Second, the extremal value of the integral of the Weyl vector is zero. Thus, no measurable 
size change occurs for classical motion, even though the Weyl vector does not have vanishing 
curl. The classical paths are precisely the ones along which no physical dilatation is ever 
measured. 

It is also interesting to note that there is a 6-dim metric, of an unexpected form that is 
consistent with collisions. It follows from the solution to the structure equations that the 
line element is of the form 

ds'^ = dx^dx^ + dx^di/i 

Therefore, if we assume that the distance ds between two particles must vanish (or nearly 
so) in order for two particles to collide, we see we must have dx"^ = 0, regardless of their 
separation in momentum, dyi. This would not be the case if we had simply imposed a 
Euclidean metric on the space. 

These results provide a satisfying unification of classical mechanics. In addition, the rel- 
ativistic version of biconformal gauging also turns out to be interesting. We have shown 
( 3|>0i) ^^^^ the method provides the best way to understand conformally invariant grav- 
ity. The results are consistent with general relativity, and improve previous conformal gravity 
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theories. The fact that we can satisfactorily gauge classical mechanics - and get something 
new - gives us a better understanding of, and more confidence in, the relativistic theory. 

There is also a suggestion of something deeper. Notice that quantum mechanics requires 
both position and momentum variables to make sense, while biconformal gauging of New- 
ton's second law gives us a space which automatically has both sets of variables. Is it possible 
that quantum physics takes on a particularly simple form in biconformal space? Anderson 
and Wheeler claim it does, deriving a path integral formulation of quantum mechanics di- 
rectly from a biconformal measurement theory 56|. It becomes possible to claim that the 
physical manifold is really a six- (or, relativistically, an eight-) dimensional place, in which 
quantum mechanics is a natural description of phenomena. 

This interpretation of biconformal space works correctly. In particular, when we use the 
covering group of the conformal group, the Weyl vector is necessarily complex. The presence 
of an "i" in the Weyl vector makes an initially real, probabilistic evolution law into a unitary 
evolution. In addition, the requirement of the scale-invariant theory for taking ratios of 
lengths to produce a meaningful measurement leads directly to the use of the product of 
probability amplitudes in computing physically measurable probabilities. Naturally, the 
proportionality between the inverse-length ^j-coordinates and momenta is taken to be 

%i = Pi 

Note that this factor drops out of Hamilton's equations, making Planck's constant classically 
unmeasurable. 

Thus, we have another way to think of quantum phenomena in a classical context. In 
this formulation, however, we have the added advantage of a direct connection to general 
relativity. It becomes possible to ask questions about quantum measurement of curved 
spacetimes in a classical context. 

VII. SPIN, STATISTICS, AND PSEUDOMECHANICS IN CLASSICAL 

PHYSICS 

A. Spin 



Now that we have a gauge theory of mechanics, we can ask further about the represen- 
tation of the gauge symmetry. A representation of a group is the vector space on which 
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the group acts. The largest class of objects on which our symmetry acts will be the class 
determining the covering group. This achieves the fullest realization of our symmetry. For 
example, while the Euclidean group ISO (3) leads us to the usual formulation of Lagrangian 
mechanics, we can ask if we might not achieve something new by gauging the covering group, 
ISpin{3) = ISU (2) . This extension, which places spinors in the context of classical physics, 
depends only on symmetry, and therefore is completely independent of quantization. 

There are numerous advantages to the spinorial extension of classical physics. After 
Cartan's discovery of spinors as linear representations of 



■ orthogonal groups in 913 (|6l|,|6^) 

M), the use of spinors for other areas of 
). Penrose developed spinor notation 



and Dirac's use of spinors in the Dirac equation 
relativistic physics was pioneered by Penrose ( 
for general relativity that is one of the most powerful tools of the field. For example, the 
use of spinors greatly simplifies Petrov's classification of spacetimes (compare Petrov 



and Penrose [65[,j68|), and tremendously shortens the proof of the positive mass theorem 
(compare Schoen and Yau ( 0| , [z^l , jzil ) Witten [z^)- Penrose also introduced the idea 
and techniques of twistor spaces. While Dirac spinors are representations of the Lorentz 
symmetry of Minkowski space, twistors are the spinors associated with larger conformal 
symmetry of compactified Minkowski space. Their overlap with string theory as twistor 
strings is an extremely active area of current research in quantum field theory (see [tsI and 
references thereto). In nonrelativistic classical physics, the use of Clifford algebras (which, 
though they do not provide a spinor representation in themselves, underlie the definition of 
the spin groups) has been advocated by Hestenes in the "geometric algebra" program j^. 

It is straightforward to include spinors in a classical theory. We provide a simple example. 
For the rotation subgroup of the Euclidean group, we can let the group act on complex 2- 
vectors, x"? = 1)2. The resulting form of the group is SU{2). In this representation, an 
ordinary 3-vector such as the position vector is written as a traceless Hermitian matrix, 

X = x'ai 
[Xf = X' [a,r' 

where at are the Pauli matrices. It is easy to write the usual Lagrangian in terms of X : 

L = ^tr (XX) - V (X) 
where V is any scalar-valued function of X. However, we now have the additional complex 
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2-vectors, x", available. Consider a Dirac-type kinetic term 
and potential 

Notice there is no necessity to introduce fermions and the concomitant anticommutation 
relations - we regard these spinors as commuting variables. A simple action therefore takes 
the form 

S = Jdt [jtr (XX) + Xa - /^X") -V{X)- XrB'aiabx'] 
The equations of motion are then 

-i tab 9V 



QX^b 

r = -i/^x" - i>^B'aiabX' 

together with the complex conjugate of the second. The first reproduces the usual equation 
of motion for the position vector. Assuming a constant vector 5*, we can easily solve the 
second. Setting x — '0^"*''*, '4> must satisfy 

This describes steady rotation of the spinor, 

The important thing to note here is that, while the spinors rotate with a single factor 
of e*^ '^, a vector such as X rotates as a matrix and therefore requires two factors of the 
rotation 

X' = e-'^^-^Xe^""-" 

This illustrates the 2 : 1 ratio of rotation angle characteristic of spin 1/2. The new degrees of 
freedom therefore describe classical spin and wc sec that spin is best thought of as a result of 
the symmetries of classical physics, rather than as a necessarily quantum phenomenon. Sim- 
ilar results using the covering group of the Lorentz group introduce Dirac spinors naturally 
into relativity theory. Indeed, as noted above, 2-component spinor notation is a powerful 
tool in general relativity, where it makes such results as the Petrov classification or the 
positivity of mass transparent. 
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B. Statistics and pseudomechanics 

The use of spinors brings immediately to mind the exclusion principle and the spin- 
statistics theorem. We stressed that spin and statistics are independent. Moreover, spin, as 
described above, follows from the use of the covering group of any given orthogonal group 
and is therefore classical. For statistics, on the other hand, the situation is not so simple. 
In quantum mechanics, the difference between Bose-Einstein and Fermi-Dirac statistics is a 
consequence of the combination of anticommuting variables with the use of discrete states. 
In classical physics we do not have discrete states. However, nothing prevents us from 
introducing anticommuting variables. In its newest form, the resulting area of study is 
called pesudomechanics. 

The use of anticommuting, or Grassmann variables in classical physics actually has an 
even longer history than spin. The oldest and most ready example is the use of the wedge 
product for differential forms 

dx Ady — —dy A dx 
This gives a grading of (— )^ to all p-forms. Thus, if a; is a p-form and r] a g-form, 

CO = a;ii...jpdx*' A • • • A dx*^ 
7] = Uij^...i^dx'^ A ■ ■ ■ A dx'" 

Then their (wedge) product is even or odd depending on whether pq is even or odd: 

Nonetheless, p-forms rotate as covariant, rank-p tensors under S0{3) (or SO (n)), in viola- 
tion of the familiar spin-statistics theorem. Under SU{2) they rotate as covariant, rank-2p 
tensors, not as spinors. 

Another appearance of anticommuting variables in classical mechanics stems from the 
insights of supersymmetric field theory. Before supersymmetry, continuous symmetries in 
classical systems were characterized by Lie algebras, with each element of the Lie algebra 
generating a symmetry transformation. The Lie algebra is a vector space characterized by 
a closed commutator product and the Jacobi identity. Supersymmetries are extensions of 
the normal Lie symmetries of physical systems to include symmetry generators (Grassmann 
variables) that anticommute. Like the grading of differential forms, all transformations of 
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the graded Lie algebra are assigned a grading, or 1, that determines whether a commutator 
or commutator is appropriate, according to 

where p,q & {0,1} ■ Thus, two transformations which both have grading 1 have anticommu- 
tation relations with one another, while all other combinations satisfy commutation relations. 

Again, there is nothing intrinsically "quantum" about such generalized symmetries, so we 
can consider classical supersymmetric field theories and even supersymmetrized classical me- 
chanics. Since anticommuting fields correspond to fermions in quantum mechanics, we may 
continue to call variables fermionic when used classically, even though their statistical prop- 
erties may not be Fermi-Dirac. Perhaps more importantly, we arrive at a class of classical 
action functionals whose quantization leads directly to Pauli or Dirac spinor equations. 

Casalbuoni pioneered the development of pseudomechanics, showing that it was possible 
to formuate an ^ — > limit of a quantum system in such a way that the spinors remain 
but their magnitude is no longer quantized ( 7^, see also Freund Conversely, the 

resulting classical action leads to the Pauli-Schrodinger equation when quantized. Similarly, 
Berezin and Marinov [t^I; and Brink, Deser, Zumino, di Vecchia and Howe Q introduced 
four anticommuting variables, 6°' to write the pre-Dirac action. We display these actions 
below, after giving a simplified example. Since these approaches moved from quantum fields 
to classical equations, they already involved spinor representations. However, vector versions 
(having anticommuting variables without spinors) are possible as well. Our example below 
is of the latter type. Our development is a slight modification of that given by Freund 

To construct a simple pseudomechanical model, we introduce a superspace formulation, 
extending the usual "bosonic" 3-space coordinates Xi by three additional anticommuting 
coordinates, 9"-, 

{r,^^} = 

Consider the motion of a particle described by [xi {t) , O"- (t)], and the action functional 



S= dt 



2 2 

Notice that 6^ = for any anticommuting variable, so the linear velocity term is the best we 
can do. For the same reason, the Taylor series in O"" of the potential V (x*, 9'') terminates: 
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Since the coefficients remain functions of x\ we liave introduced four new fields into the 
problem. However, they are not all independent. If we change coordinates from ^" to some 
new anticommuting variables, setting 



^bc - B[bc] 

where C, is an anticommuting constant, the component functions in H (^*) change according 
to 

+ (^SafbE'^iBi^ + i«£5cd + V'aC»£(.cd) xVx^ 

The final term vanishes if we choose 

while no choice of B'^^ can make the second term vanish because ijja^B^^ is nilpotent while 
^EabcB^ is not. Renaming the coefficient functions, V takes the form 

V {e') =Vo + ro" + ^^SabcB'^e'e''. 

Now, without loss of generality, the action takes the form 

S = j dt {^mx'x' + -Vo- i^aO" - ^SabcB'^e^'e' 

Varying, we get two sets of equations of motion: 

ax' ax* ax* 2 ax* 
r = itl^" + ie^B^e". 

Clearly this generalizes Newton's second law. The coefficients in the first equation depend 
only on x*, so terms with different powers of 6"" must vanish separately. Therefore, B"- and V'" 
are constant and we can integrate the 6'" equation immediately. Since [Ji^""^ = e^"' satisfies 
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we see that -B^^V element of the Lie algebra of 5*0 (3). Exponentiating to get an 

element of the rotation group, the solution for 9"" is 

The solution for depends on the force, —dV/dx^, in the usual way. 

It is tempting to interpret the 9°' variables as spin degrees of freedom and -B" as the 
magnetic field. Then the solution shows that the spin precesses in the magnetic field. 
However, notice that -B^^V in 5*0(3), not the spin group SU{2). The coordinates 6°' 
therefore provide an example of fermionic, spin-1 objects. 

One of the goals of early explorations of pseudomechanics was to ask what classical 
equations lead to the Pauli and Dirac equations when quantized. Casalbuoni see 
also [tG] showed how to introduce classical, anticommuting spinors using an — > limit of 
a quantum system. Conversely, the action 



S= dt 



1 W + _ Vo (x) - (L • S) Vls - (S ■ B) 



where L is the orbital angular momentum, S = ~^£%c^ d'^, and Vls is a spin-orbit potential, 
eads to the Pauli-Schrodinger equation when quantized. Similarly, Berezin and Marinov 
TtI ]. Brink, Deser, Zumino, and di Vecchia and Howe [3| , introduced four anticommuting 

variables, 6°' to write the pre-Dirac action, 

where 

dXn „, 



dX ' 

and a is a Lagrange multiplier. The action, SDirac, is both reparameterization invariant 
and Lorentz invariant. Its variation leads to the usual relativistic mass-energy-momentum 
relation together with a constraint. When the system is quantized, imposing the constraint 
on the physical states gives the Dirac equation. 

Evidently, the quantization of these actions is also taken to include the entension to the 
relevant covering group. 
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C. Spin-statistics theorem 

Despite the evident classical independence of spin and statistics, there exists a limited 
spin-statistics theorem due to Morgan j3] • The theorem is proved from Poincare invariance, 
using extensive transcription of quantum methods into the language of Poisson brackets - 
an interesting accomplishment in itself. A brief statement of the theorem is the following: 

Theorem: Let L be a pseudoclassical, Poincare-invariant Lagrangian, built quadratically 
from the dynamical variables. If L is invariant under the combined action of charge conju- 
gation (C) and time reversal (T) then integer spin variables are even Grassmann quantities 
while odd-half-integer spin variables are odd Grassmann quantities. 

Proof relies on extending the quantum notions of charge conjugation and time reversal. 
As in quantum mechanics, charge conjugation is required to include complex conjugation. 
For fermionic variables, Morgan requires reversal of the order of Grassmann variables under 
conjugation 

This insures the reality property [r]^*)* = rj^*, but this is not a necessary condition for 
complex Grassmann numbers. For example, the conjugate of the complex 2-form 

dz A dz* 

is clearly just 

dz* A dz 

and is therefore pure imaginary. We must therefore regard the TC symmetry required by 
the proof as somewhat arbitrary. 

t -t 
T —T 

Whether this is an allowed Poincare transformation depends on the precise definition of 
the symmetry. If we define Poincare transformations as those preserving the infinitesimal 
line element, dr, then reversing proper time is not allowed. Of course, we could define 
Poincare transformations as preserving the quadratic form, dr"^ = Qapdx'^dx^ , in which case 
the transformation is allowed. 



49 



Despite its shortcomings, the proof is interesting because it identifies a set of conditions 
under which a classical pseudomechanics action obeys the spin statistics theorem. This is 
an interesting class of theories and it would be worth investigating further. Surely there is 
some set of properties which can be associated with the classical version of the theorem. 
Perhaps a fruitful approach would be to assume the theorem and derive the maximal class 
of actions satisfying it. 

There are other questions we might ask of spinorial and graded classical mechanics. A 
primary question is whether there are any actual physical systems which are well modeled 
by either spinors or graded variables. If such systems exist, are any of them supersymmet- 
ric? What symmetries are associated with spinorial and fermionic variables? Is there a 
generalization of the Noether theorem to these variables? What are the resulting conserved 
quantities? What is the supersymmetric extension of familiar problems such as the Kepler 
or harmonic oscillator? 

The statistical behavior of fermionic classical systems is not clear. Quantum mechani- 
cally, of course, Fermi-Dirac statistics follow from the limitation of discrete states to single 
occupancy. This, in turn, follows from the action of an anticommuting raising operator on 
the vacuum: 



Since classical states are not discrete, there may be no such limitation. Do anticommut- 
ing classical variables therefore satisfy Bose-Einstein statistics? If so, how do Fermi-Dirac 
quantum states become Bose-Einstein in the classical limit? 

The introduction of pseudomechanics has led to substantial formal work on supermani- 
folds and symplectic supermanifolds. See js^, jsil and references therein. 



Clearly, the field of classical mechanics has no conclusion, and we do not provide one 
here. Within each topic we have tried to provide more questions than answers. However, in 
the process of collecting these results, we have observed a few patterns. In closing, we take 
note of those. 



a 



.t|0) = |l) 
a^a"^ = 




VIII. OBSERVATIONS 
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1. New elements in classical physics work their way into the field from fundamental 
research areas, notably quantum field theory and general relativity. The former has 
contributed spinors and anticommuting numbers while the latter lends the tools of 
differential geometry to the study of symplectic manifolds. 

2. Many of the new insights have been seen only in one or two dimensions. In these cases, 
it remains an open question whether the properties even exist in higher dimensions, 
which higher dimensions, and why in those dimensions. This applies particularly to 
the study of inequivalent Lagrangians and Bohmian quantum mechanics. 

3. Comparatively httle use is made of the classical physics ArXiv. Researchers in the 
area would benefit by using this ready reference tool. 

4. Classical mechanics is now strongly infiuenced by quantum mechanics. In addition to 
Bohmian quantum mechanics, which seeks to realize quantum physics as some sort of 

classical system, there is phase space quantization, which accomplishes much the same 
thing in a different way. An additional approach is suggested by the gauge theories 
of Section 6. These programs demonstrate broad overlap between the classical and 
quantum worlds. 

APPENDIX A: ALTERNATIVE LAGRANGIANS IN HIGHER DIMENSIONS 

We show that the Lagrangian 




where a \ ^, 6j is a time-independent constant of the motion and 9y — jv is a unit vector 
in the direction of the velocity, solves one of the Euler-Lagrange equations. 




First notice that 
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We therefore have 



Xi 



d dL 
dtdx' 



dL_ 



Xi 



dt 



X V 



e 

1 da (x, ^) 



a{x,C. ) ^ a {x, v) x' ^ (^^' 



V V 



dx' 



d^ - X* 



^2 dx' dx' 
d dL dL\ 



T 



dt dx'^ dx^ J 



x\ d r Oi{x,^) x\ da{x,v) 

Xj— r- / — d^ -\- Xi 

V dt 



X V 



^2 V dt V 

^2 dx^ dx^ 



d dL _dL\ 

dt dx"^ dx^ J 



X a\x,v) 

V I Xi h X 

V 



f 1 da a (x, v) x* 
\v dt V 
^ r 1 da{x,0 ^^ , . d df .,df 
e <9x' ^^""'dtdx' ^ dx^ 



•■ -i .i 

XjX — + f X 



da a 



—X V 



<^,xi ^ X —— - X' — 

^2 * dt 9x* dx^ 

da . f d df df 
dt \dt dx^ dx^ 



= 

Possibly the function / may be chosen so that the remaining equations of motion are satisfied. 



APPENDIX B: ARBITRARY NUMBER OF EXTREMA IN KEPLER ORBITS 



The problem of global properties of orbits remains open - power law forces have been 
studied and found to have limited numbers of extrema, but non-monotonic force laws 
allow arbitrarily many extrema. We provide a simple example here. 

Consider the potential 

V = a{r — rof^ 
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The potential V has energy and effective potential, 

'-'^ ' ^'^ + a (r - ro) 



E = -mf^ + ^ ^ ^, _ ^-^2p 

2 2mr"' 



2mr2 



+ a (r - ro 



,2p 



This effective potential has an arbitrarily strong minimum near tq. The exact location of 
the minimum is given by 



= V^ff = 2pa (r - ro) 



2p-l 



mr 



3' 



We solve this approximately as follows. Let r = ro + a. Then 



= 2par^ 1 + 3— + 3^ + „ 
u » „ r^ rg 



a 

To 



a 



m 



Now suppose a << ro so that we can neglect the ^ terms. Then in order to have solutions 
we must have 



M2 



or 



2p-l 



2pamrQ 
M2 



2B+2- « 1 (Bl) 

This may be satisfied by choosing p sufficiently large. Thus, r = ro + a is the approximate 
position of the extremum. This solution for r is a minimum since 



K// = ^ + 2p (2p - 1) a (r - ro)^^"^ > 



Now, setting r = ro + a + £, and expanding the effective potential to second order about the 
minimum at ro + a. 



eff 



2mrl ( 1 



2mrl 



. 2a ^ 
1 + 2e + 



a 



+ aa^P 1 



M2 



2pe 



+ 2p{2p-l) 



2mrl 



2a , 2n 
1 + — + ^ 1 + aa^P 
ro r^ 



+ 



2mrl 



e' + 2p {2p - 1) aa^P-h^ 
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The first term is just an overall constant, while the linear term vanishes because tq + a is a 
minimum. The third term is a harmonic oscillator potential. The approximate equation of 
motion for the oscillator is 



+ 



m \ 2mrQ 



+ 2p (2p - 1) aa2P"2 ) £ 







By eq. fjBl|) . the squared frequency becomes 
M2 



2m^rQ 



2p , 
+ - 2p 

m 



1 ) aa 



2p-2 



2m^ 







m 

M2 2p , , 

^ 2 2 + - 1) "7^ 3 

Im^TQ ma 2pamrQ 

M2 



2m'^rQ 



aro 



(2p-l) 



The frequency may be made arbitrarily large, at fixed angular momentum M, by increasing 
p. This means we may have arbitrarily many extrema per orbit. 



APPENDIX C: THE KORTEWEG-DE VRIES EQUATION 

While we have so far stayed within particle mechanics, many classical field theories also 
have interesting properties. Of particular interest are the "integrable systems" such as the 
KdV and the sine-Gordon equations. These differential equations turn out to have infinitely 
many constants of motion. 

The KdV equation is a one dimensional, third order field equation that provides a good 
example of hidden symmetries. Here we briefly examine some of its properties. The inter- 
esting history of the equation stretches over more than a century jsO]. We will begin with a 
modern form of the equation. 

Consider any function of the form u = f{z) = f {x — vt) . Substituting, we find that / must 
satisfy 

= [-Qfd + vd + d^) f 

Let g = f + c, this becomes 



= {-6gd +{v + 6c)d + d^) g 
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so choosing c = —v /6 we have simply 

= -Qggx + gx 

Integrating twice we find the quadrature, 



df 

X — vt 



These solutions for u propagate with unchanging shape / and constant velocity v. It can 
be shown that pairs of solitary waves can pass through one another and emerge unchanged. 
It has been suggested that the infinite hierarchy of constants of the motion of this system 
is related to the existence of such soliton solutions. Showing how these constants arise will 
simultaneously illustrate some techniques of classical field theory. Our treatment follows 
that of Abraham and Marsden which is recommended for further detail. 

First, we show that the KdV equation may be described as a Hamiltonian system. For 
particle motion expressed in canonical coordinates, we can define the Hamiltonian vector 
field, Xh which is everywhere tangent to the phase space motion. Conversely, the classical 
motion of the system is along the integral curves of this vector field. Restricted to any 
solution curve, Xh is therefore given by 

/ dx^ dpi 



\c={x{t),p{t))- ( 



\dt^ dt 



fdH dH 



\ dpj ' dx^ 
-\ \c={x{t),p{t)) 
We therefore can characterize Xh everywhere by writing 

— — ^ R aX^ 

or more simply using differential forms, 

dH {v) =uj{Xh,v) 

for any vector field, v. 

The same relationship holds in classical field theory. For the KdV equation, we can define 
a symplectic form as follows: 

u {u,v) = - I dx I dy [u {y) v {x) — u [x) v 
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Here u and v are arbitrary vector fields. Now suppose the Hamiltonian is given as an integral 
over a Hamiltonian density, 

[ f[u{x)]dx 



Hamilton's equations are then involve functional derivatives. For the differential of the 
Hamiltonian, 

dH (v) = / 

J —c 

so equating to the symplectic form. 



dx-J- (x) V (x) 
ou 



dH (v) =lo{Xh,v) 
dx— (x) v{x)^ - J dx J dy [Xh (y) v (x) - Xh (x) v (y)] 

We seek an expression for Xh- First, write Xh in the form 

dG 



f 

J — c 



X 



H 



dx 



Then, integrating by parts and disregarding surface terms, we obtain. 



dx^ {x) v{x) = \ I dx i dy 
ou 2 



dG dG . 

—V {X) - —V (y) 



= ^ J dxv (x) J dy 
and since v is arbitrary we have 



dG_ir ^^dG 
dy 2 J dx 



dyv (y) 



Therefore, 



Now consider the KdV equation. 



Xh^ 



d_5l 
dx Su 



The time evolution of u is given by the tangent vector field Uf. We ask if we can write this 
vector field as a Hamiltonian vector field Xh for some Hamiltonian H. Equating 



ut = Xh 
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we require 

OX du ox ^ 

du 

and it is easy to see that we can take 

H ^ J dx + ^ul 

We therefore have a Hamiltonian system, and the KdV equation may be studied in terms 
of a Hamiltonian flow. 

We now show that the KdV equation possesses infinitely many constants of motion. 
Define an infinite set of Hamiltonian vector fields and Hamiltonian densities by acting re- 
peatedly on Xi and fi according to 



= [2audx + aux + hdl) 



Su 

d 5fn+l 



dx Su 

The first expression always exists, but the second is possible as long as each new X^+i is a 
Hamiltonian flow. In order for there to exist a Hamiltonian such that 

dH{v)^uiXH,v) 

we require the integrability condition 

= d^H{v) = duj{XH,v) 

Essentially, this condition reduces to the equality of mixed partial functional derivatives. 
We omit the inductive proof that shows that the condition is satisfied for all Xn, as long as 
it holds for the initial set. Letting fi = jl it follows that 



^2 = Qx (3rt^ - lixx) = ^uux - d^u 

T, T~ = (3^* - ^xx) 

ox ou 



Sf2 o_.2 

Su 



/2 = + \ul 
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so /2 is the Hamiltonian density for the KdV equation, 



Since the inductive hypothesis holds, the entire set of Hamiltonian vector fields and 
Hamiltonian densities /„ exists. 

Finally, we are ready for the proof that there exist an infinite number of constants of the 
motion of the KdV equation. Consider the higher order "Hamiltonians" given by integrating 
the /„ : 

-f^n = y /n (a^) Ax 

We compute their Poisson brackets with one another by integration by parts. 



ex 



{Hn, H^} — f2 (Xn, X, 

= \ j j [Xn (y) Xm {x) - Xn (x) Xm (?/)] 

^\ j d^ j dy 



= / dx^X^ 
ou 



— J dx^-^ {2audx + au^. + hd^ 



m—l 



5u 



^ I dx [ 2a^—ud^ . 

ou ou 



+ a-—Ux — ; h 0— — a 



Su 



Su 



Su ^ Su 



— J dx 
+ J dx ^Ux 



Su I Su 

Sfn Sfm-1 



Su Su 



^ Su Su 



Su 



(aux + 2audx + bdl) ^ 

Su 



- / dx Xn+ 1 



Su 

{Hn+1, Hm-l} 



Now iterate this relationship. First suppose n and m are either both even or both odd. 
Then without loss of generality we take m — n — 2k > 0. Setting m — n + 2k and iterating 
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k times we have: 

{Hni Hm} = {Hn, Hn+2k] 

— {Hn+k, Hn+2k-k} 

— {Hn+k, Hn+k} 

= 

where the last step follows by antisymmetry of the bracket. Now let m — n + 2k + 1. Again 
iterating k — 1 times, and then one more time, give 

{Hn, Hnj} = {Hn, Hn+2k+l} 

— {Hn+k-1, Hn+k} 
= {Hn+k, Hn+k-l} 

But the last two lines are negatives of one another, and therefore vanish. Therefore, all of 
the Hn have vanishing Poisson brackets with one another. In particular, since H2 is the 
original Hamiltonian, {H2, Hm} — for all m, and the evolution generated by H2 leaves all 
Hm constant. Since the evolution by H2 generates solutions to the KdV equation, all H n 
are constants of integration of the KdV system. 

The KdV equation has interesting quantum properties as well. It can be shown that 
the Schrodinger equation with time-dependent potential u {t) has solutions with a fixed 
energy spectrum - the time-dependence of the potential does not change the energies of the 
soltutions. The proof hinges on the Lax theorem, which states that the KdV equation is 
equivalent to the equation 

ut = [L, A] 

where 
and 

L = dl + u 

This latter operator L is just the Schrodinger operator with potential u. The proof of the 
theorem follows by direct calculation: 

utf = [idl + 6ud^ + 3u^, dl+u]f 
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This is satisfied if —u satisfies the KdV equation. The full proof of the resulting isospectral 
theorem may be found in 5^. In light of this relationship between a remarkable classical 
system and its equally striking quantum properties, one wonders whether the relationship 
between the classical and quantum mechanics may be much like the relationship between 
the real line and the complex plane. Just as real functions often display their full character 
only when analytically extended to the complex plane, many classical systems may show 
their true natures when quantized. The KdV equation provides an excellent example of this 
- solutions to the KdV equation, when used as quantum potentials, are isospectral despite 
time-dependent potentials, and there may be a profound connection between the KdV and 
Schrodinger systems. 

The existence of equations such as the KdV equation, which have infinitely many inde- 
pendent conserved quantities is remarkable in several regards. For example, Goldstein, Poole 
and Safko observe that the KdV equation provides a counterexample to the converse of the 
Noether theorem [3| • Thus, while symmetries of an action lead to conserved quantities, the 
KdV and other equations have infinitely many conserved quantities without corresponding 
symmetries. 

There is little systematic theory of these so-called "integrable systems." In fact, we lack 
even a clear definition of this concept of integrability. Still, there has considerable recent 
progress (see, for example, |5l| and references therein). 

A related question is whether such systems exist in higher dimensions. As with many 
modern results in classical mechanics, examples are limited to one or two dimensions, and 
it is unclear whether we are seeing properties of geometries or only of the real and complex 
number systems. 



APPENDIX D: CHAOS 



While the study of nonlinear and chaotic systems is beyond the scope of this review, one 
common example provides an interesting case of a higher order differential equation for a 
classical system. The onset of chaos may be visualized by studying the fixed points of the 
logistic equation, 

Xfc+i = axk {l + hxk) 



60 



for varying values of the parameters a and b (see the review articles by May j52], May and 
Oster js^, as well as May ) As the values of these parameters change, the number of 
fixed points passes through bifurcation points, leading to more and more frequent doubling 
of the their number. At a finite value of the parameters, the number of fixed points diverges 
and the behavior of the system is said to become chaotic. 

This equation may be converted into a nonlocal equation of continuous motion for a 
one dimensional system. The fixed points of the discrete equation then become periodic 
solutions for the continuous system. Replace the discrete sequence Xk with a function x (t) , 
which must satisfy 

x{t + l) = ax (t) [1 + bx (t)] 
The left hand side may be expanded in a Taylor series as 



n! ' dt"- 

n=0 



SO at any time t, the function x must satisfy 

El d X 7 2 

— — ax — abx = 
n! c/t" 

n=0 

This certainly qualifies as a higher order differential equation! 

We can find the fixed points from the continuous representation as well as from the 
discrete one. At /c*^-order fixed points of the discrete system, we require periodicity of the 
form X {t + k) = x [t) . To examine the consequences of this condition, we employ a common 
technique for periodic systems j^]. 

Suppose Xi(t) are independent solutions to x {t + k) = x (t) . Then a general solution may 
be written as a superposition of these, so Xi (t + k) must be some superposition: 

Xi {t + k) = ciijXj (t) 
j 

Periodic solutions then satisfy 

Now diagonalize aij. If yi are the new basis functions and Aj the eigenvalues, then 



y- (t + k) = Xiyi it) 



61 



implies 

where tTj is any periodic function with period k. Now, the i/i satisfy the equations 

n=0 

n=0 

Consider the long-term behavior of yi. Suppose Aj > 1 so that yi(t) diverges at late times. 
From the first equation yi must satisfy 

-aby^ = 

as t ^ OO, so the limiting value of yi is zero, which is inconsistent. Therefore, we require 
Aj < 1 so that yi is also converging to zero at late times, and must approximately satisfy 
both 

1 d^y^ 



n=0 

OO 



EUn fin 

n=0 

Since these are linear, we may write 

y^ = e"'* 

Then 

OO / Nn 



n=0 



OO / , xn 



and we need both 



n=0 



e"-' = a 



e'^"^ = Xi = a'' 
thereby determining the (asymptotic) eigenvalues. 
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